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Abstract 



We consider a generalized convolution, linking Structure Functions (SF) 
F2 for nucleons, F 2 A for a physical nucleus and f PN ' A for a nucleus, com- 
posed of point-nucleons. In order to extract F£ we employ data on F% ,A and 
the computed f PN < A . Only for Q 2 « 3.5 GeV 2 do data permit the extrac- 
tion of F 2 A (x,3.5) over a sufficiently wide x-range. Applying Mellin trans- 
forms, the above relation between SF turns into an algebraic one, which 
one solves for the Mellin transform of the unknown FQ- We present in- 
version methods leading to the desired F<£, all using a parametrization 
for C(x,Q 2 ) = F^{x, Q 2 )/F2 (x, Q 2 ). Imposing motivated constraints, the 
simplest parametrization leaves one free parameter C{x = 1,Q 2 ). For 
Q 2 = 3.5 GeV 2 its average over several targets and different methods is 
(C(l, 3.5)} = 0.54±0.03. We argue that for the investigated Q 2 , C(x -> 1, 3.5) 
is determined by the nucleon-elastic (NE) part of SF. A calculation of the 
latter comes close to the extracted value. Both are close to the SU(6) limit 
uv(x,3.5) = 2dv(x,3.5) for parton distribution functions. 

It has been recently proposed to extract the ratio C(x, Q 2 ) = F^{x, Q 2 )/ 'F^(x, Q 2 ) of the 
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n, p structure functions (SF) from data on ratios cr( 3 H) / <r( 3 He) of inclusive cross sections. An 
electron beam with E « 11-12 GeV is expected to provide the required wide, continuous 
kinematic ranges. The proposal has already elicited discussions, mainly on the stability of 
the suggested extraction techniques for C(x,Q 2 ) 

An experiment, involving the SF of the lightest isobars with minimal I- spin symmetry 
breaking is presumably best suited to obtain C(x, Q 2 ). However, even when an upgrading of 
the beam will be approved, results on the A=3 proposal are not expected before 2006. One 
therefore wonders whether existing data on F 2 may furnish similar information. Below we 
explore this possibility and present first results for C(x,Q 2 ) and F^x^Q 2 ) from inclusive 
scattering data on D,C and Fe. 

In the following we limit ourselves to inclusive scattering of unpolarized electrons from 
randomly oriented targets A. The cross section per nucleon for beam energy E reads 

d 2 a^(E;9,u)/A _ 2 rxM 2 Ffjx, Q 2 h 

dndv ~~ M M H^ + tan {9/2) FA(x,Q 2 )\ ' (1) 

where 9 and v are the scattering angle and the energy loss. Ff 2 (x, Q 2 ) are nuclear structure 
functions (SF) per nucleon, expressed in terms of the squared 4-momentum transfer Q 2 = 
q 2 — v 2 and the Bjorken variable x = Q 2 /2Mv (M is the nucleon mass) with range < x < A. 

Henceforth nuclear and nucleon SF are assumed to be related by some generalized con- 
volution F A = f * F N , for instance ||[7| 

F h A (x, Q 2 ) = jf ^ [jf p PN ' A & Q 2 ) + -jfn N ' A (z, Q 2 )F£(x/z, Q 2 )] (2a) 

A j 

£f PN *(z,Cf)FW(x/z,Cf). (2b) 



In Eq. (H) fpn' A i s ^ ne SF for a nucleus, composed of point-nucleons, where initially a p 
or a n absorbs the virtual photon; f PN ' A is their average over the number of protons and 
neutrons in the target A(Z, N). F^ N ' stands for the similarly averaged nucleon (N) SF (We 
often drop arguments of functions when there is no danger of confusion) []. 



x * In Eqs. (0) appear admixtures of nucleon SF B, which are negligible for the involved Q 2 . 
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Eq. (Q) describes partons from nucleons in a nucleus, but not those from other sources 
(virtual bosons) neither does it account for anti-screening. Both limit the use of (fj) to x > 
0.15-0.20 @. Eq. (@) is estimated to be valid for Q 2 > Q 2 C ^ 2 - 2.5 GeV 2 

We briefly mention previously suggested extraction methods for F£. Those dealt mostly 
with D data and use f PN in the Impulse Approximation (IA) [T^jT^]. It has for instance 



been emphasized that Eq. @j is a Fredholm integral equation for the unknown F% . Dis- 
cretization in x and z produces a set of linear equations with a solution, tending to the exact 
answer for Az — > [14]. The strong variation of the 'kernel' f PN apparently hampers an 



actual application of the above. We also mention iteration methods to deconvolute nuclear 



effects for a D target [12,15 



A different deconvolution has been suggested for the proposed precision data, relating 
F A IJ. For those one may define super-ratios @,[| 



t , [Ff(x)/F£(x)} 
9 {X) a[F?(x)/FZ(x)} + b ' (4bj 

with, in principle, arbitrary a and b. For given those ratios contain C(x) explicitly, and 
are from Eq. @ seen to depend implicitly on C(x/z). Iteration determines C(x). 

Finally we recall the expansion of F lyN \x/z) in Eq. (§) around the maximum z ~ 1 of 
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F A ( X ) = F 2 (A °(x) + A^(x)[xF 2 (A ° (1) (x)] + M 2 (x)[xF^ N){1 \x) + if 1(2) (i)] + .... , (5) 

with M n (Q 2 ) = f*dzf PN > A (z,Q 2 )(l - z) n . The first ^-independent term F^ N \x) ad- 
equately represents the series for x < 0.4. Insufficient knowledge of F^ ^ix)-, the Z-th 
derivative of F^ip)-) prevents the use of the expansion (|5|) beyond that range. Although 
not the region of prime interest, the low-x result will guide the extraction of F£. 



We first recall calculations of F A from Eq. 



with Fj!' n as input. For nucleons in 



vacuum there are abundant data on F^ 



and less accurate ones for Ff [|18j. Inclusive 



data do not reach the elastic point x — 1, but the used parametrizations do. With no direct 
information on F£, one frequently uses the 'primitive' assumption 

which is the first term in Eq. (f|) for the D and thus holds only for small x. The question 
of interest is how modifications of Eq. (H) for x > 0.4 influence F A (x, Q 2 ), Eq. (|]). 

The SF f PN ' A is a many-body property and requires a model for a calculation, such as 



the perturbative I A, which to lowest order contains the single-hole spectral function [19,20 



There nucleons appear off their mass shell and a prescription is needed to handle those. We 
prefer a non-perturbative version, based on a generalized Gersch-Rodriguez-Smith (GRS) 
theory, where the SF of nucleons F N are reasoned to be on their mass shell. That method 
moreover allows a computation of f PN < A beyond its lowest order |7| JTl||2T| -p3|. 



PN,A 



The above program has initially been realized for A > 12. It has been shown that / 
is nearly independent of A, and using Eq. (^) one proves the same for F A . Data indeed 
show that for x < 0.8, the ratios /i A ' A ' = [d 2 a eA /A]/[d 2 a eA ' /A'] « F A /F A ', A, A' > 12 differ 
differ from 1 by less than 2-3% |2l||2|j2l. 

The above does not hold for the lightest nuclei, in particular not for A' — > D, when 

d 2 a eA (x,Q 2 )/A F A (x,Q 2 ) 



V A (x,Q 2 



.A,D, 



x,Q< 



d 2 a eD (x,Q 2 )/2 Ff>(x,Q 2 ) 



(7) 



are the EMC ratios, with their characteristic deviations from 1 in the range x < 0.9 [25 



For A < 4 and given iViV interaction, one may presently compute with great precision 
nuclear ground states p6| , as well as f PN > A in Eq. (|2|). Using those, GRS calculations of 
inclusive cross sections have recently been completed for D [B7f and 4 He ||28||. 



We now address the inverse problem of trying to extract F^ from data on Ff' and 
the above-mentioned computed /. First one needs to obtain nuclear SF and we shall use 
JLab data for D P9] , |3T| and for C, Fe |3T] , supplemented by some older NE3 data of more 



restricted kinematics WA - We mention two approaches: 



i) i?-ratios: Ideally one performs a Rosenbluth extraction of F A from cross sections ([I]) 
for fixed x, Q 2 and varying 9 or E, which provides 

Unfortunately, the set of (x, Q 2 ) points in the JLab and NE3 data does not enable the 
above Rosenbluth extractions. Instead one frequently invokes an empirical expression such 
as R A (x, Q 2 ) ~ 0.32/ Q 2 . It assumes A, as well as x-independence and prescribes the depen- 
dence on Q 2 . The above empirical R actually contradicts theoretical predictions |TT] . 
We shall return to this point below. 

ii) Theoretical information: Using the primitive F n ' pr , computed cross sections ([[]) 
as function of 9, v generally agree with data in the deep-inelastic regions, x < 1, (Q 2 > 



A,th 
2 



Q 2 ). There we select data, for which the relative deviation a(= a(x,Q 2 )) = (F : 
pA,expy pA,exp g^jgggg | X — ck| < 0.2 and which, for given 9 appear to change smoothly 
with x. We suggest to attribute the above deviations in equal measure to both structure 
functions F A in ([!]), which enables the definition of quasi-data 

F^ d (x,Q 2 ) = a(x,Q 2 )F A > th (x,Q 2 ) (9) 

By construction, the F A ' qd produce an exact fit to cross sections. 

Use of Eq. (§) for the purpose of inversion requires input for fixed Q 2 and running x, 
whereas cross section data are for essentially one beam energy E and cover separated x bins. 
Q 2 varies mildly within those x-bins, but not when going from one bin to a neighbouring 
one. It is thus necessary to make cuts for fixed Q 2 , containing a maximal number of x-points 
with x < 1.0 — 1.1. After careful interpolation between quasi-data, the above appears only 
fulfilled for Q 2 between 3-4 GeV 2 and then only for 0.5 < x < 1.1. It is clearly impossible 
to reliably extrapolate the manipulated quasi-data into the crucial range x < 0.5. 

At this point we invoke the empirical observation that for x$ ~ 0.16 — 0.18 , F%(xo, Q 2 ) 
and -F 2 D (xo, Q 2 ) are practically constant as function of Q 2 ||17|| . We exploit this by considering 
Eq. (Q) for small x < x$. Upon substitution of the small-x part of Eq. (||) into @ and use 
of Jq 4 dzf PN ' A (z) = 1, one easily shows that for any target 
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F*' nud {0.l8) « ^(0.18) w 0.34 (10a) 
F 2 A ' nuci (0) = F 2 <A °(0) = ^(0) ^ (10b) 

The latter extension to x = is in agreement with Eq. ([5]) and will be needed below 0. 

We recall that f PN > A in Eq. (pj) relates only to nucleons in the target and as a result, 
Eqs. fllOap hold only for the nucleonic component F < A ' nucl . That component may now be 
interpolated for x < 0.55, completing knowledge of F^' nucl (x) over the entire relevant range 
x<l.l. 

After the critical remarks on the empirical form of the i?-ratio it comes as a surprise that 
in the range 0.5 < x < 1, i 71 2 Fe,,?d (x, 3.5) and the empirically extracted Frf 6 3.5) agree 
to within ±5%, and even for 1 < x < 1.5 to within ±12%! It vindicates Arrington's claim 
that even a 100% uncertainty in the phenomenological R incurs only a 5% uncertainty in 



We proceed as follows. For Fe, for which the number of available deep-inelastic data 
points is largest, we use F Pe ' qd from method ii). With insufficient C, D data for application 
of that method, we exploit Eq. (^) and find 

rpC,R 

rpC,qd _ r 2 rpFe^qd 
r 2 — -pFe^l 

F°> qd = \pF]- l F^ qd (11) 

EMC data for \i c are in the desired range 0.3 < x < 0.80 where they require some smoothing. 
For x < 0.30 we interpolated the purely nucleonic component as in ii) and added for 0.80 < 
x < 0.88 averages of Be and Al data |[25|| . 

We now introduce Mellin transforms (MT) and their inverses, which for real g(x) are 
defined as 

POO 

g(u) = I dxx u - 1 g(x) (12a) 

(12b) 



o 



g{%) = - 

7T 



poo 

\ dtRe 


~g(a + it) ' 


Jo 





2 * The above explains why for all A and Q 2 , EMC ratios fx A (x, Q 2 ) cross 1 for x w 0.18 (25 1 
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The constant a is chosen such, that g(u) = g(a + it) is free of singularities in the complex u 
plane to the right of the imaginary u axis shifted by a. 

Application of Eq. ( |12a| ) to Eq. (0), turns the generalized convolution into a linear 
relation, which can be solved for -F 2 n (w, Q 2 ). For fixed Q 2 



F?{u) 



f PN > A (u + l)Fi N) < 



u 



+ g(N(NE)) 



r. 



■n(A) 



A 



F A (u) 



u 
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. _ —pP( u ) _ 

NfPN,A( U +l) N 2 



—Q v + g n 

N . 



(13a) 
(13b) 



with Q 



N,NE 



combinations of the standard nucleon static formfactors, making up the NE 



parts of F^ ,NE (see for instance Ref. [11]). Above we denote by the A(N, Z)-dependent 



n(A) 



right-hand side of Eq. ( |13b|) . For exact input F-, it should coincide with F^u), the MT of 



F%{x). 



Next we consider u = in Eq. ( 13b ). Since the normalization of f PN > A implies 
f PN > A (l) = 1, one finds 



A,nucl 



(« = 0) 



§mu = o) 



N 



g p + g r 



(14) 



Consequently, Eq. ( |12a| ), when used in ( |10a| ) implies that none of the MT of SF are defined 
for u — 0. This has numerical consequences also in the immediate neighbourhood of u = 0. 

Unfortunately, the direct inversion Eq. dl2b|) of F(u) runs into serious numerical prob- 
lems. We therefore take recourse to indirect methods, all featuring a parametrization of the 
ratio C in 



F?(x, Q 2 ) = F 2 "(x, Q 2 ; d k ) = C(x, Q 2 ; d k )F%(x, Q 2 ) 
C(x,Q 2 ;d k ) = Y J d k (Q 2 )(l - x) k 



(15) 



k>0 



Eq. ([lOaD implies a first constraint C(0)(= F 2 n (0)/F 2 p (0)) = E^o^^^ 2 ) = 1 Next ' one 
may exploit F 2 n (x) = F^ pr (x), Eq. (g), small x < 0.35. We use only F™(0.2,3.5) = 0.75 
which for k max = 2 leaves one free parameter do = C(l). 
It is convenient to re-parametrize F^x, 3.5) as follows 
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F%{x,Q 2 ) =ar a2 ^c m (l-x) m ;2 >0.02 (16a) 

m>l 

= 0.42 ;a;<0.02 (16b) 

A small a 2 <C 1 in Eq. (16a) produces the observed rise ofF%(x) for small x, while the cut-off 



in Eq. (16b) avoids the singularity in (16a). Eq. (10a) again implies the same for F^' n (x 



In the region 0.02 < x < 0.9, F% , Eq. ( |16a| ), practically coincides with the parametrization 
in Ref. |T7[] , both agreeing with data averaged over resonances. 

We now explore some extraction methods, comparing relevant computed or extracted 
expressions A, related to F£ and parametrized forms B: 

I) The extracted j% {A \u), E Q- (PS) and the MT of F™, Eq. (||). 

II) F A,qd (x), Eq. (H), with the nuclear SF F 2 A (x), computed from (Q), now using Eq. 

(0)- 

2 

The parameters dk, are determined by minimization of variances w(d Q ) = £\ Ai — 5, 
(or of relative variances). Table I summarizes our results. The spread of the individual 
entries C(l) from I) reflects variations in the smallest and largest u, retained in the above 
sum (see remark after Eq. (|T4"D). An empty entry indicates the absence of a minimum 
in the studied intervals. We note that in those cases the slope of the variance is very 
small and that the range of C(l) is actually compatible with values from real minima. The 
results for the two methods and for the three targets produces a well-determined average 
(C(l,3.5)) = 0.54 ±0.03. 

In Fig. 1 we show the full C(x, 3.5), as well as i 71 2 p '"(x, 3.5). On the right abscissa 
are marked the values 2/3, 3/7, 1/4 for C(l), corresponding to exact SU(6) symmetry, 
dominance of S = over S = 1 di-quark coupling and the same for the ^-component S z for 
x — > 1 (see for instance Ref. [fL5] ) 0. 

We now argue that for our purposes the regarded Q 2 ~ 3.5 GeV 2 is not 'large'. First 



3 * The extracted F% allows an evaluation of the Gottfried sumrule 5c(3.5) = 
^{dx/x)[F%{x,Q 2 ) -F£(x,Q 2 )] = 0.251. This is close to the recent value 5 G (4.0) = 0.256 ±0.026 



we note that F2 ,n (x,Q 2 ) has a first inelastic threshold N + m, m being the pion mass, or 
x thr{Q 2 ) — [1 + 2Mm/Q 2 ]~ 1 , i.e. a finite x-distance from the elastic point x — 1. That 
interval shrinks with increasing Q 2 , but for the case at hand 2^(3.5) ~ 0.93, which is 
marked by a vertical line in Fig. 1. Data show that .Ff(:E, 3.5) is negligibly small beyond 

The above implies that C(x — > 1,3.5) can be ascribed to the NE part F N,NE , i.e. to 



static N form factors (cf. Eq. (|13b|) ). Disregarding the electric form factor G n E of the n, one 
finds (/ijv are the magnetic moments of the nucleons) 



limC(x,Q 2 

X— »1 



1 + 



4M 2 (1 (Q 

Q 2 



ftp 



(17) 



with 



7(Q 2 



l^pG p E {Q 2 ) . a n (Q 2 ) G^(Q 2 )//i, 



G P M (Q 2 ) ' « p (g 2 ) G P M (Q 2 )/ \i % 



(18) 



Recent data show that 7, a p , a n deviate from 1 and equal for Q 2 = 3.5, 5.0 GeV 2 : 7=0.552, 
0.349 and a n /a p w 1.2, 1.1 [0,[3§]. Eqs. (0), (|18D then yield 



C(x = 1, 3.5) w 0.61 ; C(x = 1, 5.0) w 0.56 



(19) 



The large Q 2 limit essentially depends on the same for a n /a p . For a value 1, C(x — > 1, Q 2 — > 
00) = 0.469, which shows that in the above sense C(l, 3.5) is still far from a scaling limit. 

The above value 0.61 is reasonably close to the extracted one and has also been entered 
on the abscissa in Fig.l. Either one definitely exceeds previously cited values C(l) ~ 0.42 
from D data (cf. Ref. ||15|| ). The above seems to indicate that SU(6) symmetry u v (x, 3.5) = 
2d v (x, 3.5) for the up and down quark parton distribution functions is only mildly broken. 
It is moreover in agreement with globally extracted distribution functions B 



p7| , because of nearly identical contributions from the dominant small x region. We note in 
passing that a finite outcome requires C(0, Q 2 ) = 1. For CTEQ parametrizations from global 
parton-distribution functions C(0, Q 2 ) = (1 + 5C(0)) [pS[| , leading to a diverging Sg- 
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The above reasoning rests on forms for the nucleon SF, averageKd over resonances. 
Retaining those in detail, the sharp first inelastic threshold will become blurred, but it 
seems reasonable that C(1,Q 2 ) for relatively low Q 2 remains determined by static elastic 
form factors. That point is presently under study. A similar large-Q 2 result has been found 
from quark-hadron duality considerations, pushed to the extreme for the part F^ NE ^ BO . 

Summarizing, we have described methods to extract the neutron Structure Function from 
existing inclusive scattering data on various targets. The experimental material from several 
targets provided consistent values for C(l, 3.5). It is very desirable to plan the forthcoming 
inclusive scattering experiments on D and 4 He with the available 6 GeV beam |41 |, such 



that there be coverage of continuous x < 1-range for at least one Q 2 . Confrontation of those 
data with accurate calculations, possible for those nuclei p7| , |28|j will sharpen the present 
outcome for F 2 n (x, Q 2 ),C(x, Q 2 ). 

The authors thank G. Petratos for giving detailed information on the proposed A=3 
experiment and J. Arrington for putting at our disposal D and Coulomb-corrected NE3 data 
for C and Fe. ASR profited from discussions with G. Salme, E. Pace, several experimentalists 
at JLab and in particular with W. Melnitchuk . 
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Figures. 



Fig. 1. The ratio C(x,3.5) = F£(x, 3.5)/F|(x, 3.5) for Q = 3.5 GeV 2 from data on D, 
C, Fe. The drawn line corresponds to C(l) = 0.54 and the band represents the spread in 
the result, from averaging over different targets and methods. The vertical line for x = 0.93 
marks the pion threshold x t h r {3-5). The numbers on the right abscissa are quark model and 
QCD predictions for C(l), while 0.61 is the NE limit flUf) . Also entered are and i^, 
corresponding to C. The band for the latter is hardly noticible in F£- 
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TABLES 

TABLE I. Values of C(l,3.5) from minimizing variances w(do). Results are from data on 
D,C,Fe, using extraction methods I), II). The spread in results for I correspond to varying u-intervals 
in the variances. No entry corresponds to cases for which there is no minimum within the above 
intervals. 





D 


C 


Fe 


I 


0.55 ±0.05 




0.55 ±0.03 


II 


0.55 


0.50 
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